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I review recent results obtained within chiral effective models, on the phase structure of hot quark 
matter in a strong magnetic background. After a brief introduction, I focus on the results obtained 
within two chiral models improved with the Polyakov loop. The models differ for the content 
of interactions, but both of them are tuned to reproduce Lattice QCD thermodynamics at zero 
and imaginary chemical potential. One of them takes into account an explicit Polyakov loop 
dependence of the coupling; the other one neglects this contribution, but takes into account multi- 
quark interactions. A comparison between the phase diagrams of the two models is presented. 



X 



The many faces of QCD 
November 2-5, 2010 
Gent Belgium 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Hot quark matter in strong magnetic fields 



Marco Ruggieri 



1. Introduction 



The study of the influence of external factors on the vacuum of Quantum Chromodynamics 
(QCD), is one of the most attractive topics of modern physics. Thanks to Lattice QCD simula- 
tions, it is established that at vanishing quark chemical potential, two crossovers take place in a 
narrow range of temperature; one for quark deconfinement, and another one for the (approximate) 
restoration of chiral symmetry [|], ^ ||, Q ||]. Besides, powerful analytic and semi-analytic tech- 
niques have been developed to understand the coupling between chiral symmetry restoration and 
deconfinement, see [Q] and references therein. 

An alternative approach to the physics of strong interactions, which is capable to capture 
some of the non-perturbative properties of the QCD vacuum, is the Nambu-Jona Lasinio (NJL) 
model 101, see Refs. [|] for reviews. In this model, the QCD gluon-mediated interactions are 
replaced by effective interactions among quarks, which are built in order to respect the global 



symmetries of QCD. The NJL model has been improved recently [^, [10|], in order to be capable 
to compute quantities which are related to the confinement-deconfinement transition of QCD. In 
pure glue theory, color confinement can be characterized by the expectation value of the Polyakov 
loop [|n|], which is an order parameter for the center symmetry of the color gauge group. In 
theories with dynamical fermions, the Polyakov loop is still a good indicator for the confinement- 
deconfinement transition, as suggested by Lattice data at zero chemical potential [jl], g |3[ gl). 
Motivated by this property, the Polyakov extended Nambu-Jona Lasinio model (P-NJL model) has 
been introduced |l^, in which the concept of statistical confinement replaces that of the true 
confinement of QCD, and an effective potential describing interaction among the chiral condensate 
and the Polyakov loop is achieved by the coupling of quarks to a background temporal gluon field, 
then integrating over quark fields in the functional integral representation of the partition function. 



Chiral models with the Polyakov loop have been studied extensively in many contexts ||12|, |13|, [IJ, 

[1, mini Sim, mi. 



In this talk, I report on results obtained in Refs. [ |22| , |23i [24[ ], about deconfinement and chi 



ral symmetry restoration of hot quark matter in a strong magnetic background. Related Lattice 



studies can be found in [25, 2q]; previous studies of QCD and of QCD-like theories in a magnetic 



background can be found in [ |27| , [28| , [29[ |30| ]. Strong magnetic fields might be produced in non- 
central heavy ion collisions [31, 32]. For example, at the center-of-mass energy reachable at LHC, 
^snn ~ 4.5 TeV, the magnetic field can be as large as ' eB « 15m^ according to []32|]. It has been 
argued that in these conditions, the sphaleron transitions of finite temperature QCD, give rise to 
Chiral Magnetic Effect [31, 33, |^]. This adds a phenomenological interest to the study presented 



here, beside the purely theoretical one aforementioned. 

The plan of the talk is as follow: after a brief introduction to the model used in calculations, I 
discuss the results about Polyakov loop expectation value and chiral condensate in a magnetic back- 
ground. Then, I collect these results in the fomi of a phase diagram in the temperature-magnetic 
field strength plane. Finally, I discuss briefly about improved results that take into account a tree- 
level coupling among the chiral condensate and the Polyakov loop, and show within the model how 
this coupling modifies the response of hot quark matter to a magnetic background. 



'Here I measure eB in units of the vacuum squared pion mass ra^; then, eB = corresponds tofi~2x lO'^ Tesla. 
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2. The model: propagator and one-loop thermodynamic potential 

I consider a model in which quark matter is described by the following lagrangian density: 



^ = iif{ifD^-m)\if + ^,; 



(2.1) 



here i//^ is the quark Dirac spinor in the fundamental representation of the flavor SU (2) and the 
color group. A sum over color and flavor is understood. The covariant derivative embeds the QED 
coupling of the quarks with the external magnetic field, as well as the QCD coupling with the 
background gluon field which is related to the Polyakov loop, see below. Furthermore, 



(2.2) 



Here T correspond to the Pauli matrices in flavor space. The 4-quark term corresponds to the 
original NJL interaction [^. On the other hand, the 8-quark term was introduced in [ p2| | to stabilize 
the NJL vacuum. The interaction lagrangian is invariant under SU {2)v (E> SU {2)a U{l)v- In 
the chiral limit, this is the symmetry group of the action as well, if no magnetic field is applied. 
However, this group is broken explicitly toU(l)y®U{l)\0U{l)v if the magnetic field is coupled 
to the quarks, because of the different electric charge of u and d quarks. Here, the superscript 3 
in the V and A groups denotes the transformations generated by T3, T375 respectively. Therefore, 
the chiral group in presence of a magnetic field is U{1)y ®U{1)\. This group is then explicitly 
broken by the quark mass term to ?7(l)y. I will present results obtained within a model with G 



taken as constant and / 0: this model is named P-NJLg in [17], and it is tuned in order to 



reproduce Lattice data at zero and imaginary chemical potential. Then I will compare these results 
with those of a slightly different model, in which = but G depends on the Polyakov loop. The 



latter model is named EPNJL in []17[], and it is tuned in order to reproduce the same Lattice data 
used to tune the P-NJLg model. Therefore, both of them offer an effective description of the QCD 
thermodynamics at zero and imaginary chemical potential. It is thus interesting to compare their 
predictions on the response of hot quark matter to a magnetic background. 

To compute a temperature for the deconfinement crossover, I use the expectation value of the 
Polyakov loop, L, defined as 

L = ^Tr.exp (ipXaAl) , (2.3) 
where A4 is a static, homogeneous and Euclidean background temporal gluon field, coupled mini- 



mally to the quarks via the QCD covariant derivative [|10||. In the Polyakov gauge, which is conve- 
nient for this study, Aq = iXs<p + iX^(p^; moreover, since this work is done at zero quark chemical 



potential, L = V , which implies A\ = 0. This choice is also motivated by the study of [27], where it 
is shown that the paramagnetic contribution of the quarks to the thermodynamic potential induces 
the breaking of the Z3 symmetry, favoring the configurations with a real- valued Polyakov loop . 

I work in the one-loop approximation and neglect the pseudoscalar condensates; moreover, I 
make the assumption that condensation takes place only in the flavor channels To and T3. The mean 



field interaction term Eq. (2.2) can be cast in the form 



= -2^/1^ [GI + 2G8l^] - GI? - 3G8r^ 



(2.4) 
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where Z = —{uu + dd). The quark propagator is easily expanded in terms of Landau levels using 
the Leung-Ritus-Wang formalism [^], namely 



where Ep{x) corresponds to the eigenf unction of a charged fermion in magnetic field, and Ep{x) = 
7o(£'p(x)) ' 7o. In the above equation P = {po + iA4,0,J2y^2k\QfeB\,p,), where k = 0,1,2,... 
labels the k^^ Landau level, and J2 = sign{Qf), with Qf denoting the charge of the flavor /; A^, is 
a projector in Dirac space which keeps into account the degeneracy of the Landau levels, see [38] 
for details. 

The trace of the /-quark propagator is minus the chiral condensate (//), with f = u,d. Taking 
the trace in coordinate and internal space, it is easy to show that the following equation holds: 

(//) = -Nj^i^J ^^^(L,L,r|p„^) . (2.6) 



*:=0 

Here, 



'i^{L,Lj\p,,k) = UA-2^{L,L,T\p„k) , (2.7) 
and ^ denotes the statistically confining Fermi distribution function, 

1 + 2Lel^ "'f + Le^l^ '^f 

'^{L,L,T\p„k) = ^ , (2.8) 

where 

(oj = pl + 2\QfeB\k + M^ , M = mq + 2GL + AGsL^ (2.9) 

The coefficient ^^ = 2 — 5ko takes into account the degeneracy of the Landau levels. The function 
Ua is a smooth UV-regulator, which is more suitable, from the numerical point of view, in the 
present model calculation with respect to the hard-cutoff which is used in analogous calculations 
without magnetic field, see [23] for details. 

In the one-loop approximation, the thermodynamic potential H is given by 

a = '^{L,L,T) + Um- Z r°°^^{L,L,T\p„k) . (2.10) 

f=u,d k 

In the above equation Um = GI? + SG^L^. The last addendum arises from integration over the 
fermion fluctuations in the functional integral representation of the partition function; its contribu- 
tion is written in terms of the kernel 

^(L,L, T\p„k) = NMp)(Of + I log^ , (2.1 1) 

with 

^{L,L,T\p„k) = \ + 3Le-^'"f + 3le~^^'"f + e~^^'^f . (2.12) 
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Chital Condensate 



Thin Polyakov Loop 





Figure 1: Left panel: Cubic root of the chiral condensate, as a function of temperature and magnetic field. 
Right panel: Expectation value of the Polyakov loop, L, as a function of temperature and magnetic field 
strength. Results correspond to the P-NJLg model. 



The potential temi ^[L,L,r] in Eq. ( 2.10 ) is built by hand in order to reproduce the pure 
gluonic lattice data [|l^, [131 ], 



Analytic form of the coefficients a{T), b(T) can be found in [13]. Here it is enough to notice that 
they depend on one parameter, Tq, which corresponds to the deconfinement scale. In the pure glue 
theory, Tq = 270 MeV. Dynamical fermions induce a dependence of this parameter on the number 



of active flavors [19]. For the case of two light flavors to which I am interested here, Tq = 111 MeV. 



3. Deconfinement and chiral symmetry restoration 



In this Section, I firstly show the results obtained in [ |230 within the model with interaction 



specified in Eq. (2.2). The numerical values of the parameters can be found in the original reference. 



The model is named P-NJLg, following the nomenclature of [|17[]. Numerical data for the chiral 
condensate and the Polyakov loop expectation value are shown in Fig. |l]. They are obtained from 
the minimization procedure of the one-loop thermodynamic potential in Eq. ( 2.10| ). 

Data on the chiral condensate, left panel of Fig. |l|, are in agreement with the magnetic catalysis 
scenario: both the numerical value of the chiral condensate and the critical temperature, Tjf^, are 
increased by the magnetic field. The latter is defined as the value which maximizes \d'L/dT\. On the 
other hand, deconfinement crossover is poorly affected by the magnetic field. The deconfinement 
temperature, Ti, is identified with the temperature at which dL/dT is maximum. At = the two 
crossovers take place simultaneously at T = 175 MeV. On the other hand, at eB = 19m| the two 
critical temperatures are split of an amount of w 15%. This split is also measured within a different 
chiral model calculation [p7|], and persists even in the chiral limit [^]. 

I briefly compare the results of the P-NJLg model, with those obtained within the EPNJL 
model [ JTtI ] in which Gg = 0, but an explicit dependence of the NIL coupling on the Polyakov loop 
is taken into account (I take L = L^): 

G = g{l-aiL^-2a2L^) . (3.1) 
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Figure 2: Left panel: Chiral condensate E, measured in units of the chiral condensate at zero temperature 
and zero field, Y-q, as a function of temperature, for several values of the magnetic field strength. Right panel: 
Expectation value of the Polyakov loop, L, as a function of temperature, for several values of the magnetic 
field strength. In the figures, the magnetic fields are measured in units of m^. Results correspond to the 
EPNJL model. 



It has been shown in [39] that this kind of dependence naturally arises in the low energy limit of 



QCD, see also [40]. The numerical values of a\ and a2 have been fixed in [17] by a best fit of 



the available Lattice data at zero and imaginary chemical potential of Ref. [37], ai = a2 = o; = 
0.2 lb 0.05. The numerical values of the other parameters are standard and can be found in [24]. 

In Fig. ^ I plot data for the chiral condensate £, measured in units of the condensate at van- 
ishing temperature and magnetic field, that is Zq = 2 x (—253 MeV)^, and the expectation value of 
the Polyakov loop as a function of temperature, computed for several values of the magnetic field 
strength. The most striking feature of the results is that Ti and are tied together also in a strong 
magnetic field. At S = the critical temperatures are Tp = Tj^ = 185.5 MeV. Besides, at eB = \9m\, 
= 199 MeV and = 201 MeV. This is due to the fact that the NIL coupling constant in the 
pseudo-critical region in this model decreases of the amount of the 15% as a consequence of the 
deconfinement crossover. Therefore, the strength of the interaction responsible for the spontaneous 
chiral symmetry breaking is strongly affected by the deconfinement, with the obvious consequence 
that the numerical value of the chiral condensate drops down and the chiral crossover takes place. 
The picture remains qualitatively and quantitatively unchanged at eB = 30nr^. In this case, Tl = 224 
MeV and = 225 MeV. 

4. The phase diagrams 

Data on deconfinement and chiral symmetry restoration pseudo-critical temperatures are col- 
lected in Fig. Ill, in the form of a phase diagram in the eB — T plane. The left panel corresponds 



to the EPNJL model [)24|]; the right panel corresponds to the results of the P-NJLg model []23|]. In 
the figure, the magnetic field is measured in units of m|; temperature is measured in units of the 
deconfinement pseudo-critical temperature at zero magnetic field, namely Tb=o = 185.5 MeV for 
the EPNJL model, and Tb=o = 175 MeV for the P-NJLg model. In the Figure, the golden region 
corresponds to a phase in which hot quark matter is deconfined, but chiral symmetry is still broken 
in a non-perturbative way by the chiral condensate. 
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Figure 3: Left panel: Phase diagram in the eB—T plane for the EPNJL model. Temperatures on the 
vertical axis are measured in units of the pseudo-critical temperature for deconfinement at eB = 0, namely 
Tc — 185.5 MeV. Right panel: Phase diagram in the eB—T plane for the PNJLg model. Temperatures on the 
vertical axis are measured in units of the pseudo-critical temperature for deconfinement at eB = 0, namely 
Tc = 175 MeV. In both the phase diagrams, T^, Ti correspond to the chiral and deconfinement pseudo- 
critical temperatures, respectively. The golden region denotes the portion of phase diagram in which hot 
quark matter is deconfined and chiral symmetry is still broken spontaneously. 



The most astonishing feature of the phase diagram of the P-NJLg model is the entity of the 
split among the deconfinement and the chiral restoration crossover. The difference with the result 
of the EPNJL model is that in the former, the entanglement with the Polyakov loop is neglected in 
the NJL coupling constant. However, the two models share one important feature: both deconfine- 
ment and chiral symmetry restoration temperature are increased by a strong magnetic field. They 
disagree quantitatively, in the sense that the amount of split in the EPNJL model is more modest 
in comparison with that obtained in the P-NJLg model. A similar phase diagram calculation has 
been performed within the Quark-Meson model in [27|, see also the talk by Ana Mizher in this 
conference. In [27] it is shown that the magnetic field leads to the split of deconfinement and chiral 
symmetry restoration crossovers, if vacuum fermion fluctuations are properly taken into account; 
furthermore, both and are enhanced by a magnetic field, in excellent agreement with the 
behavior of the NJL-like models. 



The most recent Lattice data about this kind of study are those of Ref. []25p. In [pSQ, the 
largest value of magnetic field considered is eB 0.75 GeV^, which corresponds to eB/m\ 38. 
The Lattice data seem to point towards the phase diagram of the EPNJL model, since no split is 
measured among the two crossovers. On the other hand, the results of [25] might not be definitive: 
the lattice size might be enlarged, the lattice spacing could be taken smaller (in [ p5| ] the lattice 
spacing is a = 0.3 fm), and the pion mass could be lowered to its physical value in the vacuum. 
As a consequence, it will be interesting to compare the model results with more refined data in the 
future. 



5. Conclusions and outlook 



In this talk, I have reported on recent results [ 22 , 23 , 24 ] about the phase diagram of hot quark 
matter in a strong magnetic background. I have considered two models: the P-NJLg model and 
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the EPNJL model, following the nomenclature of [17], which differ for the interactions content, 
but both are tuned to fit Lattice data about QCD thermodynamics at zero and imaginary chemical 
potential. I have also compared the results with those of [ 27 1 in which a Quark-Meson model is 
used. The most striking similarity among the several models, when the contribution of the fermion 
determinant is properly kept into account, is that they all support the scenario in which chiral 
symmetry restoration and deconfinement temperatures are enhanced by a strong magnetic field. 
The models differ quantitatively for the amount of split measured (very few percent for the EPNJL 
model, and of the order of 10% for the other models for eB ^ 20m|). Furthermore, I have also 
compared the model results with those obtained on the Lattice [25|. 

It is fair to admit that the study presented here might have a weak point, namely, it misses a mi- 
croscopic computation of the dependence of the NJL coupling on the Polyakov loop. Equation (3. 1) 
is just a particular choice of such a functional dependence. Different functional forms respecting 
the C and the extended Z3 symmetry are certainly possible, and without a rigorous derivation of 
Eq. (3J_), it merits to study the topics discussed here using different choices for G{L) in the near 
future. For these reasons, I prefer to adopt a conservative point of view: it is interesting that the 
EPNJL model, which is adjusted in order to reproduce Lattice data at zero and imaginary chemical 
potential, predicts that the two QCD transitions are entangled in a strong magnetic background; 
however, this conclusion might not be definitive, since there exist other model calculations which 
share a common basis with EPNJL, and which show a more pronounced split of the QCD tran- 
sitions in a strong magnetic background. More refined Lattice data will certainly help to discern 
which of the two scenarios is the most favorable. 

The results collected here suggest several continuation lines. Firstly, it is worth to compute the 
magnetic susceptibility of the chiral condensate [41] within chiral models, which allow to treat self- 
consistently the spontaneous chiral symmetry breaking in a magnetic field. Besides, the technical 
machinery used here can be easily applied to a microscopic study of the spectral properties of 
mesons in strong magnetic field, in particular of p mesons following [p9|]. Furthermore, it would 
be interesting to make a complete analytical study of the chiral limit, to estimate the effect of the 
magnetic field on the universality class of two-flavor QCD. Finally, the phase diagram with chiral 
chemical potential deserves further study, following [E2I |43|]. 
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